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To clarify the relation between chiral symmetry breaking and color confinement, we investigate 
the Polyakov loop in terms of the Dirac eigenmodes in SU(3) lattice QCD. We analyze the low- 
lying (IR) and UV Dirac-mode contribution to the Polyakov loop, respectively, using the Dirac- 
mode expansion method. In the confined phase, the Polyakov loop (Lp) remains almost zero 
and Z3 center symmetry is thus unbroken, even after removing low-lying Dirac-modes, which 
are responsible to chiral symmetry breaking. In the confined phase, the Polyakov loop (Lp) also 
remains almost zero by UV Dirac-modes cut. In addition to the confined phase, we analyze the 
Polyakov loop in the deconfined phase and its temperature dependence. The behavior of the 
Polyakov loop (Lp) is found to be almost unchanged by the cut of low-lying or UV Dirac-modes 
in both confined and deconfined phases. 
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1. Introduction 

Nowadays, quantum chromodynamics (QCD) is established as the fundamental theory of the 
strong interaction. However, the non-perturbative properties of QCD are not fully understood, es- 
pecially, on color confinement and chiral symmetry breaking. It is interesting issue to investigate 
the correspondence between these non-perturbative phenomena [[l[ ^, ^, ^ ||, ^, ^ ||]. The lattice 
QCD calculations show the simultaneous chiral and deconfined phase transitions at finite tempera- 
ture [^, which suggests a close relation between confinement and chiral symmetry breaking. 

As for chiral symmetry breaking, the chiral condensate {qq) is directly connected to the Dirac 
operator in QCD. The chiral condensate is proportional to the Dirac zero-mode density as 

(^^) = -lim lim7r(p(0)), (1.1) 
with the Dirac spectral density p(A), which is known as the Banks-Casher relation [|l^]. The Dirac 



zero-modes are also related to the topological charge via the Atiyah-Singer index theorem [11]. 

Therefore, it is interesting to investigate color confinement in terms of the Dirac-operator prop- 
erties. Using Gattringer's formula the Polyakov loop was analyzed by the Dirac spectrum sum 
with twisted boundary condition on lattice ^, |6|]. In our previous studies |8|], we developed 
the Dirac-mode expansion method for the link-variable, and analyzed the role of the Dirac mode 
to the Wilson loop and the interquark potential. As for the hadron spectra, it is reported that the 
hadrons still exist as the bound state even without chiral symmetry breaking by removing low-lying 



Dirac-modes [12, 13]. 

In this paper, based on the Dirac-mode expansion method [^, ^], we investigate the role of the 
Dirac mode to the Polyakov loop in both confined and deconfined phases at finite temperature in 
SU(3) lattice QCD. In Sec. 2, we briefly review the formalism of the Dirac-mode expansion method 
in lattice QCD. In Sec.3, we show the lattice QCD results. Section 4 is devoted for the summary. 



2. Formalism 

In this section, we briefly review the Dirac-mode expansion method in lattice QCD [^, ^, and 
formulation of the Dirac-mode projected Polyakov loop. 

2.1 Dirac-mode expansion in lattice QCD 

In lattice QCD, the Dirac operator E> = y^D^ is expressed as 

1 4 

using the link-variable U^{x) € SU(A'^c) and the lattice spacing a. Here, we use the convenient 
notation of U-^ {x) = Uji{x — jX), and fx denotes for the unit vector on lattice in -direction. In this 
paper, we adopt hermitian y-matrices, i.e., 7^ = 7^, in the Euclidean space-time. Thus, the Dirac 
operator ID is anti-hermitian, and the Dirac eigenvalues are pure imaginary number. We introduce 
the normalized Dirac eigenstate \n) which satisfies 

ID\n) = iXn\n), (2.2) 
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with the eigenvalue A„ G R. The Dirac eigenfunction i//,, (x) defined by 

Wn{x) = {x\^) (2.3) 

satisfies £)i//„ = iXnYn- 

We consider the operator formalism in lattice QCD |^ |8||. The link- variable operator is 
defined by the matrix element 

{x\U^\y) = U^{x)d,+f,^y. (2.4) 
The Dirac-mode matrix element {n\U^\m) is expressed as 

{n\U^\m) = ^{n\x) {x\U + fi) {x + fi\m) 

X 

= Y^'^l{x)Uii{x)\ifm{x + fl), (2.5) 

using the link- variable U^{x) and the Dirac eigenfunction y^fnix). 

Using the completeness relation L„ |«) («| = 1, any operator O can be expanded in terms of the 
Dirac-mode basis \n) as 

d = ££|?i)(«IO|m)(m|, (2.6) 

n m 

which is the mathematical basis of the Dirac-mode expansion [Q, ^. 

Based on the expansion in Eq.(2^), we introduce the Dirac-mode projection operator P as 

P=Y,\n){nl (2.7) 

with the Dirac eigenstate \n), and arbitrary set of eigenmode subspace £/. For example, the IR and 
the UV Dirac-mode cut are given by 

•^«|>AiR |A„|<Auv 

respectively, with the IR/UV cutoff parameter. Air and Auy. We define the Dirac-mode projected 
link-variable operator as 

Ui;=PU^P= £ £ \n){n\U^\m){m\, (2.9) 

with the projection operator P. Using the projected link-variable f/^, we can analyze the individual 
contribution of each Dirac eigenmode to the various quantities, such as the Wilson loop ||]. 

2.2 Polyakov loop operator and Dirac-mode projection 

Next, we formulate the Dirac-mode projection of the Polyakov loop. Hereafter, we consider 
the periodic SU(3) lattice of the space-time volume V = x Nt with lattice spacing a. In lattice 
QCD operator formalism, the Polyakov-loop operator is defined by 

1 ^' 1 

^^-^n^4=^t/f (2.10) 
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with the temporal link-variable operator tj^. Taking the functional trace "Tr", we obtain the stan- 
dard Polyakov loop (Lp) as 

Tr Lp = i-Tr{n(74} = T^ti-£(x,^inf/4lx,0 

(=1 X,t !=1 

= —irY^{xj\Uj,\x,t + a){x,t + a\lJ^\x,t + 2a)---{x,t + {Nt - l)a\Ui^\x,t) 

= ^tvY^UA{x,t)U^{x,t + a)---UA{x,t + {Nt-\)a) = {Lp), (2.11) 

where "tr" denotes the trace over SU(3) color index. 

Using the projection operator P, we define the Dirac-mode projected Polyakov loop {iFp"^') as 

(LP) = ±Tr{YlU^} = l^Tr {PU4PU4P---PU,P} 

= 3^^ L («l|^4|«2)(«2|f/4|«3)--- («W,|f/4|«l)- (2.12) 

In particular, we consider the IR/UV Dirac-mode projected-Polyakov loop as 

(Lp)iR = -i-tr 52 {ni\U4\n2) ■ ■ ■ {nN,\U4\ni) , (2.13) 
|A„,|>A,R 

(Lp)uv = z^tr ^ {ni\U4\n2)---{nN,\U4\ni), (2.14) 

|A„,|<Auv 

with the IR/UV eigenvalue cutoff, Air and Auv- 

3. Lattice QCD calculation 

In this section, we study the Polyakov loop in terms of the Dirac-mode in SU(3) lattice QCD at 
the quenched level. We use the LAPACK package for the full diagonalization of the Dirac operator 
[jl4|]. We use the Kogut-Susskind (KS) formalism for reduction of the computational costs [0, ||]. 

3.1 The confined phase 

First, we analyze the Polyakov loop properties in the confined phase. Here, we use 6^ lattice 
with j8 = 5.6, which corresponds to lattice spacing a ~ 0.25fm ||]. The total number of KS 
Dirac-modes is X Aff X 3 = 3888. Fi gure |l| shows the lattice QCD result for the Dirac spectral 
density p (A ) and IR/UV-cut Dirac spectral density, 

Pir(A)=p(A)0(|A|-Air), Puv(A) = p(A)0(Auv-|A|), (3.1) 

with Air = 0.5a^^ and Auv = 2.0a^^. Both mode-cuts correspond to removing about 400 modes 
from full eigenmodes. 

We show in Figs]^a)-(c) the scatter plot of the original Polyakov loop (Lp), (Lp)iR for low- 
lying Dirac-mode cut with Air = 0.5a^^ and (Lp)uv for UV Dirac-mode cut with Auv = 2.0a^\ 
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Figure 1: The Dirac spectral density in SU(3) lattice QCD on 6^ at j3=5.6, i.e., a ~0.25fm. (a) The original 
spectral density p(A). (b) Pir(A) for IR-cut with Ajr = O.Sa^'. (c) P\i\{X) for UV-cut with Auv = 2.0fl"^ 
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Figure 2: The scatter plot of the Polyakov loop in the confined phase in SU(3) lattice QCD on 6^ at J3=5.6, 
i.e., a ~0.25fm and T = \ /{Nta) ~ 0.13GeV. (a) The original Polyakov loop (Lp). (b) (Lp)ir for low-lying 
Dirac-mode cut with Air = 0.5fl^^. (c) (Lp)uv for UV Dirac-mode cut with Auv = 2.0a^^. 



respectively. As shown in Fig.^a), the Polyakov loop satisfies (Lp) ~ 0, which indicates the con- 
fined phase. 

By removing low-lying Dirac-modes, chiral symmetry breaking is effectively restored ||8, 10, 
1^, [T^. Actually, this IR Dirac-mode cut of Air = 0.5a"^ ~ 0.4GeV corresponds to about 98% 
reduction of the quark condensate around the physical region ntq ~ 5MeV [js]]. However, as shown 
in Fig.^b), the Polyakov loop (Lp)iR remains almost zero, which means unbroken Z3 center sym- 
metry. This result indicates that the single-quark energy is still extremely large, and the system 
remains in the confined phase even without chiral symmetry breaking. In the UV Dirac-mode cut, 
the chiral condensate is almost unchanged, and the Polyakov loop (Lp)uv also remains almost zero, 
as shown in Fig.^c). These results in Figs.^(b) and (c) show that the Polyakov loop is insensitive 
to the IR/UV Dirac-mode cut. 



3.2 The deconfined phase at high temperature 

Next, we study the role of the Dirac mode in the deconfined phase at high temperature. Here, 
we use 6^' x 4 lattice at j8=6.0. The total number of KS Dirac eigenmodes is L^xN,x 3=2592. 

Figure ^ shows the original Polyakov loop {Lp), (Lp)iR for low-lying Dirac-mode cut with 
Air = 0.5a^^ and (Lp)uv for UV Dirac-mode cut with Auv = 2.0a^^, respectively. These mode- 
cuts correspond to removing about 200 modes from full eigeimiodes. As shown in Fig.||(a), the 
Polyakov loop has a non-zero expectation value (Lp) ^ 0, which shows the center group Z3 struc- 
ture on the complex plane. This property indicates the deconfined phase. 

After removing low-lying or UV Dirac-modes, as shown in Figs.||(b) and (c), the Dirac-mode 
projected Polyakov loop (Z./))ir/uv still shows the non-zero value and the Z3 center structure, which 
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Figure 3: The scatter plot of the Polyakov loop in the deconfined phase in SU(3) lattice QCD on 6^* x 4 
at j3=6.0, i.e., a ~0.10fm and T = 1 /{N,a) ^ O.SGeV. (a) The original Polyakov loop {Lp). (b) (L/.)ir for 
low-lying Dirac-mode cut with Ajr ~ 0.5a^'. (c) {Lp)\jy for UV Dirac-mode cut with Ajjy = l.Oa^K 
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Figure 4: The j3-dependence of the Polyakov loop | (ip)iR/uv I in SU(3) lattice QCD on 6^ x 4. (a) | (Lp)ir | 
for IR Dirac-mode cut with Air = 0.5fl^' and l.Oa^'. (b) |(Lp)[jv| for UV Dirac-mode cut with Auv = 
2.0fl^' and l.la^^. The original Polyakov loop |(Lp)| without cut is added. 



indicates the deconfined and Z3 broken phase. In fact, in both cut cases of IR and UV Dirac modes, 
no drastic change occurs on the Polyakov loop, apart from a constant normalization factor. The 
Dirac-mode seems to be insensitive also for deconfinement properties of the Polyakov loop. 

3.3 j3 -dependence of Dirac-mode projected Polyakov loop 

We also investigate the j3 -dependence of the absolute value of the Dirac-mode projected 
Polyakov loop, |(ip)iR/uv|» fixed A/, and L. Here, we use 6^ x 4 lattice with j8 = 5.4 ~ 6.0. 

Figure ^ (a) shows |(Lp)ir| with Air = 0.5a^^ and I.Oa^^ and Fig.^ (b) shows |(Lp)uv| with 
Auv = 2.0a^' and IJa^^. In terms of the removed number of Dirac modes. Air = 0.5a^^ and 
l.Oa^^ approximately correspond to Auv = 2.0a^^ and IJa^^, respectively. We have also added 
the original (no Dirac-mode cut) Polyakov loop [(Lp)|, which shows the phase transition around 
j8 = 5.6 ~ 5.7. Both IR and UV Dirac-mode projected Polyakov loop | (Z.p)ir/uvI show the similar 
j(3 -dependence as the original Polyakov loop \{Lp)\, apart from a constant normalization factor. 



4. Summary and concluding remarks 

In this paper, we have analyzed the direct relation between the Dirac eigenmodes and the 
Polyakov loop in SU(3) lattice QCD calculation at the quenched level. Using the Dirac-mode 
expansion method, we have carefully removed the relevant ingredient of chiral symmetry breaking 
from the Polyakov loop. 
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In the confined phase, we have found that the Polyakov loop remains almost zero even without 
low-lying Dirac-mode. These low-lying modes are relevant for chiral symmetry breaking, as the 
Banks-Casher relation indicates. However, the Polyakov loop does not show any drastic changes, 
which indicates the system still remains in the confined phase. This result is consistent with the 
Wilson loop analysis, which shows the area law and the linear interquark potential even after re- 
moving low-lying Dirac-modes ||]. We have also checked the UV Dirac-mode contribution to 
the Polyakov loop. By removing UV Dirac modes, the Polyakov loop remains almost zero. Thus, 
there seem to be no specific Dirac modes essential for the Polyakov loop in the confined phase. 

In addition to the confined phase, we have also analyzed the Polyakov loop properties in the 
deconfined phase at high temperature. In the deconfined phase, the Polyakov loop has a non-zero 
expectation value, which distributes in Z3 direction in the complex plane. Even by removing low- 
lying or UV Dirac-modes, the behavior of the Polyakov loop (ip)iR/uv seems almost unchanged, 
apart from a constant normalization factor. 

These lattice QCD results suggest no direct connection between chiral symmetry breaking and 
color confinement through the Dirac eigenmodes, which indicates that one-to-one correspondence 
would not hold between them in QCD. If it is the case, the QCD phase diagram would exhibit more 
richer structure by mismatch of chiral and deconfinement phase transitions. 
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